Charm production by cosmic muons 
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Narrow muon bundles in underground detectors permit to study muoproduction 
reactions that take place in the surrounding rock. We analyze the relevance of a 
QED+QCD reaction, muoproduction of "open charm". The contribution to double 
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muon events is estimated to be 4 — 8% of the one due to QED "trident" process, for an 
ideal detector located under a rock depth of 3 km water equivalent, and an observation 
Q4' threshold of 1 GeV. 

In recent years, there has been a certain experimental [1], |2| and theoretical interest 
^ H> [S H) H 011 "narrow muon bundles" (multiple muons with a lateral separation less than 

a few meters) in underground detectors. These events have been observed as a peak at 
small lateral separation, and interpreted as an induced flux. In fact, the energetic muons 
that propagate underground in roughly ~ 1% of the cases interact and produce other muons. 
Thence, an analysis of these events requires to consider high energy muoproduction processes, 
in the rock surrounding the detector. 

Up to now, the process considered was the muon "trident" reaction |7|] /i Z — > fi Z yU + /i~, 
where a muon pair is formed in the field of the nucleus^. For muons propagating in high Z 
materials an amplification factor Z 2 /A (= 5.5 for standard rock, A=22 and Z—ll) is present, 
due to coherent character of the reaction. This interpretation has been pursued since the 
first evidences obtained in cosmic ray experiments 0. The trident reaction leads mostly 
to narrow bundles of three or two muons in an underground detector (one produced muon 
may stop before reaching the detector); a reference ratio in an ideal detector is of 3 double 
*e-mail: vissani@lngs.infn.it 

1 It is assumed that an effective rejection of muoproduced n^, 7, e^ ... can be achieved. 



muons per triple muon, for a threshold of observation Eh = 1 GeV, and a depth h = 3 km 
w.e. of standard rock. Recent studies P, [|, however, suggest that existing interpretations 
are insufficient to quantitatively account for the whole "narrow muon" data set. 

In this work we analyze the role of another high energy process as source of prompt 
muons: production of charmed states due to cosmic (atmospheric) muons, whose relevant 
energies range fromP] E ~ 100 GeV up to tenths of TeV's (for studies in laboratory, see 
PH|). More specifically, we are concerned with the "open charm" reaction of muoproduction: 
fiN — > ficcX (X denotes a byproduct which does not concern us). This process is stipulated 
by QED and QCD interactions, while weak interactions provide the instability of charmed 
states: c — > X c — > \iX. 

In order to obtain a simple estimate of the flux of double muons due to this process, we 
adopted the collinear approximation, considering only how the initial muon with energy E 
branches into those of the final muons (E' and E") and proceeded in the following way: 

1) First, we calculated the cross section of muoproduction da^^^ccxl 'dE'dE c at leading 
order (LO) in a s , double differential in the energies of the scattered muon, E' and of the 
charm, E c (see appendix). This can be done with a limited effort by following the calcula- 



tions documented in where the cross section integrated over the hadronic phase space 
d&hadr was obtained: In fact, neglecting the gluon mass, the differential expression is simply 
d^hadr = dE c j '(8ivE^), where _E 7 = E — E' is the energy of the virtual photon emitted by 
the muonP]. In the actual calculation, that requires integrating over the photon virtuality 
Q 2 and the gluon momentum fraction x, we use the GRV98 gluon distribution [fT2] , and set: 
m c = 1.5, 1.35 or 1.2 GeV. We multiplied the cross section by the factor K(E) = (Tnlo/o'lo 
(where a is the total cross section) to describe next-to leading order QCD effects |TB|, [TJ], |T5J. 
The differential cross section increases with E' with a u l/v 2 behavior" and than rapidly 
decreases to zero in the range of energies of interest; instead, it is rather mildly distributed 
in E c . The total cross section a increases with E, due to the smaller values of x that are 
probed by the virtual photon, and to well known characteristics of the gluon distribution. 
Its value is 4 x 10~ 32 cm 2 when E — 1 TeV for m c = 1.35 GeV (almost equal to the trident 
cross section at the same energy); LO cross section increases by 50% if m c is lowered to 1.2 
GeV, while decreases by 30% if m c is 1.5 GeV. 



2 We neglect the energy loss in the rock of the charmed hadrons X c , for a 200 GeV D ± meson travels on 

average just 3 cm in the rock before decaying. 

3 Also, we found convenient to relate E c to the zenith angle and velocity of emission in the gluon-gamma 



c.m.s. as follows: E c /E 1 = (1 + (3* cos#*)/2, where (3* = y/T— 4to;?/ (p + q) 2 (p and q are the momenta of 
the gluon and of the virtual photon) 



2) We estimate a "scaling" probability dP^^/dw = BR C ^^ x p c _^{w) that a charm 
yields a muon with a certain energy fraction w = E"/E c , by first hadronizing the charm into 
a D meson (using the normalized distribution of with ed = 0.135) and then letting it 
decay with a distribution (that is, retaining only the D mass, and neglecting the Q 2 
dependence of the form factors). The resulting normalized probability p c ^n(w) falls strongly 
with the energy fraction w; the median of the distribution is in fact (E") = 0.15 x E c . We 
took as an effective branching ratio of charm into muons the value BR^^ = 8% |17[| , and 
multiplied the result by two, to account for the fact that a charm or an anticharm can yield 
a muonF]. Notice, incidentally, that the corresponding yield of triple muon is negligible, due 
to an a priori 4% suppression factor. 

3) At this stage, we can calculate the cross section da^^^j dE' dE" , where E" is the 
energy of the produced muon, and, with that, the cross section^ ct^n^^^E, f) for production 
of two muons, each with a fraction of the initial muon energy greater than /. Due to the 
behaviors of da^-^nccX an d dP c ^^ with E' and E" mentioned above, this cross section 
diminishes dramatically with /; when E — 1 TeV, it drops down by one order of magnitude 
already when / pa 0.07. This cross section enters the elementary yield of double muons in 
the detector, which depends linearly on the infinitesimal depth crossed dh! (in gr/cm 2 ): 

dY^^(E, ti) = dh' x N A x a^^^E, f) where / = ^ (1) 

Na = 6.023 x 10 23 is the number of nucleons in 1 mole (multiplying by dh', we obtain 
the density of targets per cm 2 ). The energy losses are evaluated in continuous energy loss 
approximation, Ey = (Eh+e) exp[ (h — h!) / h ]—e, where e ~ 600 GeV and ho ~ 2.5 km w.e. 
are phenomenological parameters, and Eh = 1 GeV is the (typical) threshold for detection. 
Multiplying this by the single muon flux differential in dE, dF^, we get the differential double 
muon flux induced by "open charm" reaction at the depth h. The total flux is then: 

rh roo 

F^(h)= / dF„{E,ti)xdY^{E,ti) (2) 

where we integrated over the depth of production h', and the energy E of the primary muon 
at this depth (namely, where the reaction takes place); E was related to the energy at 
the surface E in the continuous energy loss approximation, which permits us to evaluate 
dF^(E, h')/dE by using the (approximate) expression for the flux at the surface given in ITS 



4 Existing underground detectors do not distinguish between "same charge" and "opposite charge" double 
muon events. 

5 We consider only those events whose vertex is not contained in the detector. Those events profit of a 
large effective target mass, and correspond, in a sense, to the celebrated neutrino-induced single muon signal. 



(namely, dF^/dEo = 0.14 x E 2 ' 7 x ... /(cm 2 s sr GeV)). The results are shown in the figure, 
for muons arriving from the vertical direction. 
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Figure 1: Flux of narrow double muons due to open charm formation (thick curves, for 
m c = 1.5, 1.35, 1.2 GeV from lower to higher one) and to the trident reaction (thin curve, 
calculated with the cross section given in 

The contribution of open charm reaction is not large; for instance, at a depth of 3 km 
w.e. it is just 4 — 8% of the one due to the trident process. Equivalently, it can be compared 
with the flux of single muon: we get F^/F^ = 0.7, 1, 1.4 x 10 -5 in the case of open charm 
reaction (for m c = 1.5,1.35,1.2 GeV) while F^/F^ = 1.8 x 10~ 4 in the case of trident 
reaction. For an ideal detector, it would require to accumulate several hundredth (trident 
narrow double muon) events, to become statistically interesting. The smallness of the result 
has to be attributed to the relatively small value of BR^^, and to the effective leakage of 
energy of the virtual photon, during the conversion 7* — > c — > D — > [i (while for tridents, 
7* immediately materializes into muons). However, this contribution is not negligible if one 
aims at reaching the precision of 5 — 10% in the predictions. 

The following remarks illustrate other aspects of this result: 

(a) going to shallower depths, the double muon flux increases, though less rapidly than the 
single muon flux: in fact, the effective target increases with the depth (but, of course, also 
the background increases); 

(b) conversely, in deeper sites the relative contribution of the open charm process becomes 
more important, due to more energetic primaries-.^ increases (but there are practical limi- 
tations, due to the time of data taking and area of the installation); 

(c) keeping the depth fixed, and changing the angle of observation, there is an increase of 



moving toward the horizon, directly related to the increase of (but the actual geometry 
of the rock in the underground site has an essential role in practical considerations); 
(d) in water or ice, the trident curve would reduce by ~ 1.5 in comparison with the open 
charm one, due to the Z 2 /A factor. This would somehow emphasize the open charm contri- 
bution (but it should be reminded that no plan exists to have an underwater detector, with 
large area and capable to achieve a good discrimination at small lateral distances). 

In conclusion, it seems to us rather difficult to account for a large fraction of narrow 
muon bundles on the basis of the open charm process. Thus, the chances of studying heavy 
quark physics with existing underground (or underwater) detectors are quite limited. This 
result, however, adds motivations for further search of unexpected sources of background 
(and, possibly, new sources of prompt muons) when we recall the difficulties to understand 
existing data on narrow muon bundles. For the future experimental perspectives, we consider 
interesting the possibility to achieve energy and charge identification of the muons in un- 
derground detectors, as a possible handle to separate various components in a narrow muon 
bundles data set (for instance, we found that the average energy of the parent muon-that 
forms two muons through the open charm process-is rather large, above 1 TeV). However, 
even if it will be possible to obtain sufficient statistics and control of the systematics, an 
attempt to proceed further and extract a signal of production of heavy quarks from studies 
of narrow muon bundles will need more refined calculations: To accurately describe the NLO 
effects [14 , hadronization, and decay of charmed states [ 19 , |2(| ; but also those effect in the 



muon propagation, that are necessary to model the lateral distribution of the events in the 



underground detectors |2T|, |22], gj. In fact, the relatively large transverse momenta p± ~ m c 
that result from charm production and decay lead to larger lateral separations than those 
due to the trident reaction, and this could be of interest to characterize the charm induced 
events. 
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Appendix: Formulae for LO cross sections 

The LO cross section, differential in y = E^/E and z = E c /E^ is: 

da„ N ^- cX _ o? J dhgQ2 J dhgx a ^ 2) Q{x ^ 2 



dydz 9ME . 

(3) 

2m l , A df T ( 2m\ . \ df^ 

1 " # + 2 g(y) ) if + ( 1 ~ # + 6 9(a) ) Tz 



where M = 0.938 GeV, m M = 0.106 GeV and a = 1/137; a s is the strong coupling, G 
the gluon distribution function, and fi 2 the factorization scale (/i 2 = 4m 2 . + Q 2 in present 
calculation). In order to describe nuclear effects in "standard rock" nuclei, we followed [J24I , 
and weighted the gluon distribution function by the density: r s r (x) ~ 1.26 x x 073 x (1 — 
0.3 x); see eq. 1 and fig. 2 in that workQ. The adimensional functions introduced in eq. ^| are 

g(y) = (i - y)/y 2 and 

^ = 9i-9o + %c9i ~ f {92 + 2 gi ) - 2x Q ( gi - fg 2 ) + 2x 2 Q g x 

(4) 

^ = 2x Q (go-f gi )-2x 2 Q go 



where g n = 1/ z n + 1/(1 — z) n for n = 0, 1, 2; a; c = 2m 2 /pq and = Q / (2pg) are bound by 
Xmin to be lower than unity. The non-trivial limits are: (z) x mi „ = (m 2 g\{z) + Q 2 )/ (2ME 1 ), 
which results from the kinematics of the 7*g — > cc process (setting the gluon mass to zero), 
and (ii) Q 2 min ~ m 2 /g(y) which results from setting to zero the scattering angle in the 
laboratory frame; the limits on y and z, and Q 2 max follow by consistence. Notice that the 
cross section can be easily integrated analytically over z, which amounts to replace: 

OR* 91 81* its 

9o^2f3, ^_21og— -, 92 -> 1 _ (/3 , )2 (5) 

and, also, #1(2) — > pi(l/2) = 4 in the expression of x min ; the resulting expression for 
da^N^^x/dy is equivalent to the one shown in the appendix of fllfl . The cross section 



that enters the expression of the double muon yield (eq. |lj) is obtained as: 

(TnN^w(EJ) = 2 x BR C ^ x f f dy [ dz da ^^ ccX x f dw p c ^(w) (6) 

Jf Jf/v dvdz Jf/(vz) 



The last integral corresponds to an integral probability, and can be tabulated separately to 
simplify the calculation. 



6 Notice that here x is the momentum fraction of the gluon (the particle that feels the nuclear effects), as 
contrasted with xp = Q 2 / (2Pq), the variable introduced to parameterize the structure function F cc (xp, Q 2 )- 
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